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The Pomeron contribution to p p and p p scattering in AdS/QCD 



We consider the difFerential and total cross sections for proton-proton and proton-antiproton 
scattering in the Regge regime from the point of view of string dual models of QCD. We argue that 
the form factor which appears in the differential cross section is related to the matrix element of the 
stress tensor between proton states and give a procedure for computing the strength of the coupling 
of the Pomeron trajectory to the proton. We compute this coupling in the Sakai-Sugimoto model 
and find excellent agreement with the data at large s and small t. The form factor can be estimated 
in the Skyrme model or in AdS/QCD models and gives a stiffer form factor than the commonly 
used electromagnetic form factor, in agreement with our fits to data. Our model is also in good 
agreement with the measured ratio of real to imaginary parts of the forward scattering amplitude 
at large s. 

PACS numbers: 11.25.Tq, 11.15.Pg, 11.55.Jy, 13.85.Lg, 



In recent years, holographic QCD (hQCD) models based on the AdS/CFT correspondence have developed as a 
useful framework for understanding the structure of strongly coupled QCD. They have proved remarkably successful in 
reproducing meson masses, decay constants, and couplings. The difficulty of performing string calculations on curved 
backgrounds has largely limited their application to low-energy processes which can be studied in the supergravity limit 
of the string dual. However, such interesting (and experimentally relevant) processes as the high-energy scattering 
of hadrons lie outside this limit. In particular, at large center of mass energy ^/ s and small momentum transfer t, 
perturbative QCD fails - but considering the exchange of only the lowest energy confined states in the supergravity 
limit of hQCD is also insufhcient. Instead, one should include the full tower of higher spin string excitations having 
the appropriate charge and parity. In this paper, we propose a technique inspired by Regge theory and the structure 
of string scattering amplitudes to model the high s, low t limit of proton-proton and proton-anti-proton scattering. 

The outline of this paper is as follows. In section II we review the connection between Regge theory and string theory 
and discuss the AdS/QCD interpretation of the Pomeron. In section III we use the structure of QCD string duals to 
develop a model for p — p or p — p elastic scattering in the Regge limit. This model depends on four parameters: the 
slope and intercept of the Pomeron trajectory, the strength of the coupling of the Pomeron to the proton, and a mass 
scale which determines the relevant proton form factor. In section IV we compute the mass scale, the Pomeron-proton 
coupling, and the mass of the lowest particle on the trajectory, in the Sakai-Sugimoto model. In section V we fit our 
model to data and compare the fit values of the parameters with the computed values. In section VI we conclude and 
discuss some directions for future research. 



In this section we discuss some of the connections between Regge theory and string theory. While all of the results 
are known, we hope that this quick summary proves useful to string theorists unfamiliar with Regge theory and to 
experts on Regge theory who may not have kept up with the latest developments in string theory and AdS/QCD. 
It will also serve to establish some of our key assumptions and to compare them with some of the phenomenological 
literature on Pomeron exchange. While Regge theory can certainly be constructed independently of string theory, 
many of its results are more easily understood from a string-theoretic point of view; the Regge behavior of hadronic 
processes also hints strongly at the existence of a string dual description of QCD. As we will see, string duals of QCD 
illuminate some aspects of Regge theory which in the past were determined purely phenomenologically. 

We are concerned in particular with the forward behavior of proton-proton or proton-antiproton scattering. In the 
Regge limit of small t and large s, the total cross section, the difFerential cross section, and the p parameter are given 
in terms of the invariant scattering amplitude {A^P{s,t) or A^^{s,t)) ^62j as 
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I. INTRODUCTION 



II. REGGE THEORY, STRING THEORY, AND THE POMERON 
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QCD is strongly coupled in this regime, so a perturbative expansion in the QCD coupling cannot be used to compute 
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these quantities. It is natural, then, to try to compute them in the l/Nc expansion, which often provides insight 
into the properties of QCD when standard perturbation theory fails [TJ[2]- At large Nc, the QCD spectrum includes 
infinite towers of narrow mesons and glueballs with arbitrarily high spin. For example, it is easy to construct gauge 
invariant operators with spin J which have the same quantum numbers as the lo meson. These operators create higher 
mass and spin versions of the lo, which also have nonzero couplings to baryons. Since the effective coupling between 
mesons is gcs ^ l/^/N^, these states become arbitrarily narrow as Nc — > oo. Tree- level i-channcl exchange of such 
a spin J meson implies an amplitude which behaves as s"^ . Naively, terms with larger J become more important at 
large s; the sum over all such exchanges would then lead to amplitudes which grow rapidly at large s, violating the 
requirements of unitarity. 

Regge theory instead sums up these exchanges of particles with fixed quantum numbers and higher mass and spin 
(Regge trajectories) into a form that can be associated with a t-dependent pole at J = a{t) in the complex angular 
momentum plane [3]. There is a great deal of evidence that these Regge trajectories are well-approximated by a linear 
function, a{t) = a{0) + a't. At positive t (corresponding to the crossed channel for a 2 ^ 2 scattering process) this 
is evidenced by the fact that mesons can indeed be arranged into groups with J = a{0) + a'M"^ . For example, the 
trajectory which includes the / = 1 mesons p,a2, P3,a4^ has a{0) ^ 0.53 and a' ~ 0.88 GeV~^. At small negative t, 
one can extract a(t) from the differential cross section, which in Regge theory has the characteristic form 

Data from a variety of scattering processes shows that the linearity of trajectories extends to moderate values of 
negative t with the same slope and intercept as at positive t. While there is evidence that more complicated singularities 
are also required to describe the full range of hadronic data (e.g. Regge cuts) the basic picture of exchange of single 
Regge poles accounts for the structure of many hadronic scattering processes in the Regge limit. 

High-energy total cross section data quickly makes it clear, however, that a theory based purely on the exchange 
of known meson states is not sufficient to describe the scattering of hadrons. The Regge behavior described above 
leads to total cross sections atot oc (lm^(s,0))/s which behave as s"(°'~^ at large s. Since all the known meson 
Regge trajectories have intercepts q;(0) < 0.6, they cannot explain the experimental fact that total cross sections for 
hadronic processes tend to grow very slowly with increasing s, at a rate independent of the quantum numbers of the 
scattered particles. This was apparent even in the early 1960's, and led [8l[9] to propose the existence of a new Regge 
trajectory dubbed the "Pomeron," with even signature, vacuum quantum numbers, and intercept a-p{0) ~ 1, that 
governed the large s behavior of total cross sections. There have been many attempts over the years to fit total cross 
sections with a combination of Reggeon and Pomeron contributions. See for example [1111 fTTj . both of which conclude 
that the Pomeron has an intercept a-p{0) ~ 1.06 — 1.08 and a slope ~ 0.25 GeV~ . 

However, the existence of the Pomeron trajectory and its structure has remained elusive for a number of reasons. 
First of all, it is not easy to identify the corresponding trajectory at positive t, though many candidate states with the 
correct quantum numbers do exist. It is not even clear that the Pomeron trajectory is unique. One could certainly 
imagine that there are several Pomeron trajectories, just as there are a variety of Reggeon trajectories. Second, the 
gO.08 gj-owth of total cross sections implied by fits to the Pomeron intercept eventually violates the Froissart-Martin 
bound, which requires that total cross sections grow no faster than log^ s as s —^ oo. Thus single Pomeron exchange 
cannot be the full story for sufficiently high s. Whether we have already reached a value of s where such effects must 
be included is quite controversial. Attempts to test the question of single Pomeron exchange versus unitarized models 
of multi-Pomeron exchange (or extrapolations of perturbative QCD using fits to forward data) have unfortunately led 
to inconclusive results [T2]. Finally, at increasing \t\ one encounters features in the data which are not easily described 
by single Pomeron exchange with the same trajectory as at small This has led to the idea of a "hard Pomeron" 
which in some papers is treated as a separate trajectory and in others is viewed as a change in the behavior of a single 
trajectory as \t\ increases. 

String dual models of QCD [HI [H US [TS] which grew out of the AdS/CFT correspondence [HIIIHIIII] have shed 
some light on these issues [20]. In hQCD, meson states are dual to open strings and glueball states to closed strings. 
Analysis of the glueball spectrum in string duals [5TJ as well as lattice gauge theory studies [2Y support the idea that 
there is a single Pomeron trajectory consisting of glueball states, with the lowest mass state on the leading trajectory 
being a 2++ state and with the lowest mass 0++ state lying on a "daughter trajectory" (a trajectory having the 
same quantum numbers as the "leading" trajectory, but a smaller intercept). The analysis of [20] also establishes a 
connection between the soft Pomeron and the hard BFKL Pomeron, which emerges from perturbative QCD [23,. For 
theories like QCD, which exhibit confinement and logarithmic running couplings, it leads to a theory where a discrete 
spectrum of poles at positive t evolves continuously into a set of closely spaced poles on a fiatter trajectory (see e.g. 
Fig. 11 in (20j). The transition between these two behaviors is expected to occur at scales of order Aqcd- This 
picture very closely resembles one which emerges from models which generically use a fifth dimension, r, to encode 
the energy scale of the theory, and thus the running of the string tension [25' . The full amplitude therefore becomes a 
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sum over densely spaced trajectories with a{0) and a' depending on the energy scale r, so that a different trajectory 
dominates for at each value of t. Analyzing the Veneziano and Virasoro-Shapiro amplitudes in this context (see e.g. 
[251I17|) reveals a similar dependence of the trajectory slopes to that found in [5U]. Unfortunately, the region of most 
interest for our analysis, —-^qcd < ^ ^ 0; appears to be extremely model-dependent and difficult to analyze in all 
such models. Based on the picture outlined above, we will assume that there is a single Pomeron trajectory and 
focus entirely on the regime in which this trajectory promises to dominate. We also assume that one of the primary 
effects of the curved space background in AdS duals is to shift the slope and the intercept of the leading closed string 
Regge trajectory from their flat space values. Computing this shift in the region of most phenomenological interest 
is a difficult problem. For various approaches see [201 HI]. We will later present an analysis of the effective trajectory 
extracted from data that suggests a linear trajectory in the region —0.6 GeV^ < t < 0, and will fit its slope and 
intercept. 

We now review how the most important elements of Regge theory indeed emerge from string amplitudes. Here 
we consider only flat space amplitudes, which are sufficient to illustrate Regge behavior. In the next section we will 
assume that amplitudes in curved space retain much of this structure but with shifted values for the parameters of 
the Regge trajectory. 

Let us begin with an amplitude for open string exchange, which should correspond to the exchange of Reggeons. 
The Veneziano amplitude, first introduced as a model for n + n ^ lo + tt scattering, and later recognized as defining 
open string scattering, can be written in terms of the amplitude 

.Von / .^ _ T[n - ao{s)]T[m - ao{t)] 
The crossing symmetric combination appearing in Veneziano's original paper ^ is 

Aj!i,y{s, t) + 4- (t, ^) + 4r:i2} (4) 

while the open bosonic string four-tachyon amplitude is simply 

A(P1,P2,P3,P4) (X Aj^%^o}i^,t) + AjS%,o}{u,t) + AjS%,o}{s,u) (5) 

where we have ignored Chan-Paton factors and an overall constant. In the above we take ao{x) to be a linear 
function: ao{x) — ao(0) + a'^x. To avoid notational confusion later on, we use ao{x) and adx) for the linear functions 
which appear in open and closed string amplitudes and reserve the notation a(0) and a' for the intercept and slope 
of the straight line relating J to Af^ on a Regge trajectory. Thus for closed string Regge trajectories we have 
J = ac(0) + a'^M'^ and we define a^x) = adO) + a'^x. This will become clearer in the following when we contrast 
the Regge limits of open and closed string amplitudes. 

In the Regge limit (s +00 and t fixed), the Veneziano amplitude does not behave smoothly: it has poles for large, 
positive, real s. This can be remedied by giving s a small imaginary part, which corresponds physically to giving the 
meson states a small width, thus moving their poles slightly off the real axis. One finds that for large, complex s 

4oAO}(^'^) - (-«»'^°<*^r[-a„(i)] = °(*)«.)'^°(*)r[-a„W], (6) 

4o.o,o}K*)-^««)°°^*^r[-a„(i)], (7) 

(using u ^ — s in the Regge limit) and that ^^"^"(u, s), which does not have any t-channel poles, vanishes exponentially 
in s in this limit. From this we can easily see that the open string amplitude explicitly reproduces the Regge form of 
the differential cross section in ([2|. 

The linear combinations of amplitudes which are even or odd under the exchange u s have Regge limits 

^ 4o.o.o}(^>i) + ^m^M (1 + e-"«°(*))«s)'^°(*)r[-a„(t)] 

-^o" = 4o.o,o}(^'*) - 4o."o,o}(".^) - (-1 + (*))«s)'^°(*)r[-a„(i)] (8) 

The prefactors in round brackets are known in Regge theory as "signature factors" and have zeroes when ao{t) is an 
odd integer in A'^ or an even integer in A'^ . Using the product formula 

7a<,(t) °o -a„(t)/n 

n-^o{t)] = — - n r-u\ — T (9) 
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we see that A'^ has poles at ao{t) = 2k with residue ^ s^^ for k a nonnegative integer. Thus corresponds to 
the exchange of a Regge trajectory of particles with even integer spin. Similarly, A~ corresponds to exchange of a 
Regge trajectory with odd integer spin. Since in cither case we find the angular momentum of the particles being 
exchanged near a t-channel pole is J = Oi{t), for the open string we can identify = a'^ and ao(0) — ao(0) and thus 
ao{t) — ao{t). When one considers the scattering of string states with higher spin, or extends the bosonic string to the 
superstring, the general structure of tree-level amplitudes in terms of ratios of Gamma functions remains the same, 
the only difference being the addition of a kinematic factor iiro(l,2,3,4) which depends on the polarization tensors 
or spinor structures of the scattered particles and their momenta. This implies that the higher mass and spin states 
being exchanged have couplings and vertex factors tightly constrained by duality in terms of the couplings and vertex 
factors of the lightest exchanged states. 

We have now seen that the open (bosonic) string amplitude exhibits the correct pole and residue structure to 
reproduce the form of the phenomenological Regge amplitude, and also naturally includes signature factors for ex- 
change of resonances with even or odd spin. Of course, the open bosonic string amplitude has obvious problems as a 
candidate for describing scattering of mesons. In particular, the ground state, corresponding to the first pole of 
is a tachyon and the first spin one state exchanged in A~ is a massless gauge field. However, one can easily modify 
the string theory amplitude for phenomenological purposes [SI IZ] by using 

instead of ^^0,0,0} ^^d by using 
the phenomenologically determined values of the slope and intercept for meson trajectories. 

Let us now turn to closed string scattering, dual to the exchange of a trajectory of glueball states (i.e. the Pomeron). 
The simplest example is the scattering in flat space of four closed string tachyons, given by the Virasoro-Shapiro 
amplitude 

,f ^ T[-aM]n-ac{t)]n-ac{u)] 

j,,(PuP2,P3,Pi) Y^^^^^t) „ a,(s)]rha,(t) - a,{u)]T[~a,{u) ~ a,{s)] ' ^'"^ 

When we use this amplitude for phenomenological purposes we will include a kinematic prefactor _R'c(l, 2, 3, 4) and 
fit the slope and intercept to data rather using the values from the critical bosonic string. 

To take the Regge limit it is useful to write the u dependence in terms of s and t. For the simplest case of 2 ^ 2 
scattering of particles with equal mass m we have 

ads) + ac{t) + ac{u) = 4a'^m^ + 3ac(0) ee x (11) 

so that 

Ac{Pl,P2,P3,Pij = 77^ , s-iT^r iT^r , s y 12 

T[-ac{t) - ac{s)]T[ac{t) - x]r[ac(s) - x] 



Then, using the limits 



^Jf^^^->(-<^)-«, (13) 
r[-ac(s) - ac{t)\ 



and 

Tjads) + acjt) - x] 
r[ac(s) - x] 

we have the Regge limit of the Virasoro-Shapiro amplitude 



«s)^=W, (14) 



Ar'{Pi,P2,P3,P,) = W(a;.)^--(^) y[ (15) 

This amplitude has poles at ac{t) = n for n = 0, 1, 2, • • • with residue ~ s^". This corresponds to t-channel exchange 
of a Regge trajectory with J = 2n = 0, 2, 4, • • • . Note that it is not necessary to do any projection or add amplitudes 
to get exchange of only even J poles; this arises as a consequence of the fact that the leading closed string Regge 
trajectory contains only even spin states. We also note that J = 2ac(t), implying that a'^ — 2a'^ and ac(0) — 2ac(0). 

The above discussion suggests two possible amplitudes to describe the exchange of an even signature Regge tra- 
jectory whose lowest state has spin Jg = 2: (1) the Regge limit of either the even signature part of the Veneziano 
amplitude for open string theory or (2) the Regge limit of the Virasoro-Shapiro amplitude for closed string theory. 
Rewritten in terms of the actual Regge trajectories ao(t) and ac{t), including the kinematic prefactor, they are 

AT^is,t) = [Ko(l,2,3,4)(a'„s)^2]g^„a„(t)/22c^g(^^^(i)/2)(a;s)"°(*)r(2-a„(i)), (16) 
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and 

Af'^Hs,t) = [K^il, 2, 3, 4) (a^./2)-^]e— W/^(a:,./2)"-(*) ^[^^j,^/^ -T- x] - ^^^^ 

By assumption the factor in square brackets approaches a constant at large s. 

The total cross section and the parameter p = KeA{s,0)/ Im A{s, 0) are identical in either case. Provided the same 
Regge trajectory is used, they differ only by an overall constant at t = 0. However, the amplitudes do predict distinct 
differential cross sections as these involve behavior at non-zero t. Regge amplitudes are often written in the form 
(s/so)"(*) and we see that the reference scale sq differs in the two amplitudes. They also differ in their t dependence 
with the most dramatic difference being the existence of zeroes in the open string amplitude at ao{t) = —1, —3, —5, • • • 
(these are often called EXD zeroes in the literature and arise from the need to project out the even spin states from an 
Exchange Degenerate trajectory with both even and odd spins) while the closed string amplitude has zeroes coming 
from the poles in the denominator at adt) = 2^ + 2, 2%, 2x — 2, • • • . 

In the application of this formalism to proton-proton scattering discussed in the following section we use the closed 
string amplitude to model Pomeron exchange, since the Pomeron trajectory is a closed string trajectory in holographic 
duals of QCD. The situation is a bit murky, however. In dual models the proton is neither an open string nor a closed 
string but rather a solitonic excitation of open strings (as expected from large Nc reasoning [28^), so while the 
exchanged states are indeed closed strings, the states being scattered are not. 



III. A GENERAL MODEL FOR p-p ELASTIC SCATTERING IN THE REGGE REGIME 

We now use the behavior of the string amplitudes sketched above to develop a method for computing proton- 
proton scattering in a holographic string dual of QCD. One should note that at present all the proposed duals have 
serious limitations; we also lack the technical tools to calculate the full tree-level string amplitudes in a curved space 
background, which should actually govern the behavior of proton-proton scattering. As a result we will have to make 
certain approximations. Our first assumption (consistent with factorization in Regge theory) is that we can split the 
calculation into two parts. (1) We determine the vertex which governs the coupling of the Pomeron to the proton, 
which we assume is dictated by the vertex for the lowest state on the Pomeron trajectory, the 2++ glueball. Using this 
vertex, we compute the amplitude for tree-level exchange of a spin 2 glueball. (2) We then convert this amplitude into 
the Regge limit of a full tree-level string amplitude by "Reggeizing" the propagator. This is a heuristic procedure, 
described in detail below, which gives an answer consistent with the general principles of Regge theory. It should be 
a good starting approximation if the main effect of curved space and background fields on the string theory is to shift 
the slope and intercept of the Regge trajectory from their flat space values. 



A. The glueball coupling 

Our first task is to determine the coupling of the glueball field to the proton. The glueball field can be treated as 
a second-rank symmetric traceless tensor ft,^^. Old ideas of tensor-meson dominance, [121 130] applied now to the 2+"'" 
glueball rather than to the /2 meson, suggest that h^^ should couple predominantly to the QCD stress tensor T'"': 

Sint d^xh^,,T^'r (18) 

Assuming this is true, the glueball-proton-proton vertex is determined by the matrix element of the stress tensor 
between proton states 

(p',s'|V(0)|p,s) . (19) 
Using symmetry and conservation of T^^ this matrix element can be written in terms of three form factors I31j as 

(p',s'|r^,(o)b,s) = s(p',s') 

where k — p' ~p' , t ^ k"^ and P ~ {p+p') /2. The fact that the proton has spin 1/2 and mass nip implies the constraints 
A{0) — 1 and -B(O) = 0. We will see in the following subsection that the contribution from C{t) is suppressed in the 
Regge limit and that the contribution from B{t) is small compared to that from A{t). We note that gravitational 
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FIG. 1: The Feynman diagram for the exchange of a glueball in the t-channeL 



form factors of nucleons have been studied in the holographic context in |32j and also play an important role in the 
analysis of deeply virtual Compton scattering |33j : for our purposes, however, a simple analysis of their behavior in 
the Regge limit will suffice. 

How are the two ingredients - the glueball and the proton - manifest in a string dual description? Any hQCD model 
necessarily involves a theory containing gravity in a 5-dimensional (5d) space (often along with additional compact 
directions), where one of the coordinates is dual to the energy scale of QCD. There is therefore inevitably a graviton, 
which gives rise to a mode transforming as the desired 2++ glueball field. In order to have dynamical mesons and 
baryons (rather than just a baryon vertex as in [34J) there must also be fields in the theory which are dual to operators 
constructed out of quark fields. In particular, there must be a gauge field dual to the axial-vector current operator 
which creates pions in QCD. In the large Nc limit of QCD baryons may be treated as Skyrmions, that is as solitons 
of the pion field. Dual models lend themselves to this interpretation of baryons as Skyrmions, though we could also 
contemplate adding in fields dual to baryon operators in QCD, as suggested in the recent work of [35 . 

Fluctuations Hm n of the 5d background metric by definition couple to the 5d stress tensor for the matter fields via 

i J d^x^g T^T'h'''' (21) 

where ^mtv'*''^ includes a contribution from the fields which are dual to the pion. In this case, Tmn is the energy 
momentum tensor on a solitonic solution representing the proton. To reduce this to a 4d coupling of the glueball we 
expand the spin 2 piece of the metric perturbation h^'^ in terms of the glueball wavefunction, and the matter fields in 
terms of the 4d pion and vector meson fields. There is no guarantee that this will yield a 4d coupling of the glueball 
predominantly to the 4d stress tensor of the proton. However, analogy to similar results for vector meson dominance 
in dual models [36j [37] lead us to believe that this is indeed the case, and we show by explicit calculation that this 
is true in the Sakai-Sugimoto model. Making this assumption, then, we can evaluate the coupling A as an overlap 
integral of pion and glueball wavefunctions. 

Given a semi-classical solution representing a baryon it is then straightforward following the discussion in [3S] to 
compute the relevant form factors. We should note that we are computing this vertex in the large Nc limit. On the 
dual string theory side, this is a classical limit with the string coupling — > 0, so the calculation can be done in 
terms of a semi-classical solution to the equations of motion. 



B. Tree level glueball exchange 

Let us now use the generic form of the vertex discussed in the previous section to calculate the cross section due to 
glueball exchange. We will then "Reggeize" the glueball propagator wherever it appears in the amplitude, to include 
in our model scattering of higher spin glueballs. 

Consider a massive, spin-2 glueball exchanged in the t-channel (the dominant channel in the Regge limit). The 
Feynman diagram for this process is shown in FIG. [l] 

We use pi and p2 as incoming momenta and and p4 as outgoing momenta, with s, t, and u defined in the usual 
way, with k = pi ^ p^ — p2 ~ Pi the momentum of the glueball. The massive spin-2 propagator (as given in |39j ) is 



where nig is the mass of the glueball and the indices contracted at one side of the propagator are a, (3, the indices at 
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the other end are 7, S, and 
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Using the glueball-proton-proton couphng in the form of eq. (20 1 from the previous section, the ampHtude becomes 



A = 



t — Trig 



A(t)(Si7"u3)(pi + psf + ^(pi + P3fkp{uia"Pus) + ^{u^u^){k^kP - rh)] (24) 

ZTOp nip J 



v4(t)(z227^^^4)(p2 +P4)' + '4^{p2+PiYk^{u2a'^'^Ui) + ^(u2Zi4)(fc'^fc' " V'''t) 
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Using the Dirac equation we can see that 



k^^iuiYuz) = ki_,{u2"f^Ui) = 



and 



k^,{Pl+PzT = k^,{p2+Ps.T = 0. 



(25) 



(26) 



Furthermore, the second structure in the vertex will vanish when dotted into pairs of fc^. This means we can ignore 
all terms in the propagator except those of the form 7777. We will be taking the Regge limit of the amplitude, so we 
can drop any terms that will be suppressed by factors of t/ s or rrip/s. The amplitude then becomes 
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The cross section can then be calculated from the amplitude, and it will be proportional to 
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(28) 



where we suppress all terms subleading in </s or rrip/s. 

The form factor B{t) is zero at t = and slowly varying. Together with the factor of ^^^2 this implies that at small 

t, the term proportional to B^{t) is very small compared to the A^{t) term. The part of the cross section proportional 
to A(t) dominates, allowing us to drop all other terms and simply associate a XA{t) to each vertex in the amplitude. 



givmg 



da _ X^s^A'^it) 



dt 7r(t - m?Y 



(29) 



as the cross section for spin 2 glueball exchange. Note that the appearance of in the numerator is precisely the 
correct dependence expected from the exchange of a spin-2 particle. The denominator {t — to^)^ will be replaced 
by the square of the "Reggeized" propagator we describe in the next subsection. 
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C. Reggeizing the propagator 



Having computed the cross section for exchanging the lightest state on the 2++ trajectory, we must include the 
higher spin states on the trajectory, which correspond to stringy excitations on the curved background. As noted 
above, the computation of the full string amplitude is prohibitively difficult. Instead, we analyze in greater detail the 
scattering amplitude for four closed strings in flat space, and from this extract the propagator for a ^-channel closed 
string exchanged in the Regge limit. 

As discussed in the previous section, the scattering amplitudes for closed bosonic strings and for closed superstrings 
take the form 

, T[-a{t)]T[-a,{u)]T[-a,{s)] 

•^^W, 77^ rTTFr 77\ TTTR r~\ fTT^^^^^i'^^^, • • ■ ) (30) 

1 [-ac(t) - ac[s)\L [-ac[t) - ac(u)\l [-ac(u) ~ ac\s)\ 

where is a kinematic factor with no poles, which depends on the momenta and polarizations of the scattered 
particles, acix) is some linear function related to the spectrum of closed strings: 

aj^x) = ac(0) + cL^x. (31) 

We assume that this basic form of the amplitude holds true in (weakly) curved space, with ac(0), and the kinematic 
factor undetermined and dependent on the details of the geometrical background. This amplitude boasts many 
features that are important to modeling Pomeron exchange. For example, it is completely symmetric under the 
exchanges of s,i, and u. We know experimentally that proton/proton scattering and proton/anti-proton scattering 
have the same behavior in the limit where Pomeron exchange dominates, and we know that crossing symmetry 
therefore requires that the amplitude be completely symmetric under exchanges of s,t, and u. In addition, the 
amplitude has the correct pole and residue structure to describe the exchange of even spin, vacuum quantum number 
states. 



To find the proper "Reggeized" replacement for the spin 2 glueball propagator in ( 29 1 , we first expand the amplitude 
around one of the i-channel poles, which occur at a^it) — n, with rt = 0, 1, 2, . . . : 

{l + ac{s))---{n + a^{s)){-x + a^{s))---{-X + n-l + arXs)) 
^ " n!r[n-x](n-ac(t)) ■ ■ ■ ^ ^^'^ 



where x is defined as in eq. (11) and is evaluated at ac{t) = n. As we would expect for the sum of exchanges 
of all particles on the Regge trajectory, the residue of the pole is a polynomial in s, and we can identify the leading 
behavior of this polynomial as s"' if the particle being exchanged has spin J. That is 

_ P2n(a^g)(.sV(t,e,) + ...) 

nlT[n - x]{n - a,{t)) ^^^^ 

where 

X,(pi,P2,-..)«sV(e») + --- (34) 

with /(ei) some unknown function of the polarizations of the scattered particles which will factor out of our calculation 
in the end. P2„(a^s) is a polynomial of degree 2n whose first term is (a^s)^". The spin of the nth particle on the 
trajectory is then 

J=k + 2n. (35) 

We can therefore use the leading s behavior of the kinematic factor Kc to arrange the spin of the lowest particle 
on the trajectory to be whatever we want, after which the higher spins are completely determined. In our case, the 
n = pole should correspond to a spin-2 particle, meaning that k = 2. This implies that the trajectory of particles 
contributing to the amplitude in the Regge limit (where s is large) consists only of even spin particles, consistent with 
the Pomeron coupling identically to particles and their anti-particles. When we assume that the Pomeron is the dual 
of a closed string, this comes about quite naturally. By contrast, if we had assumed that the Pomeron is dual to an 
open string, we would have used the Veneziano amplitude, which has poles for both even and odd spin particles. 

Let us now now relate our string amplitude parameters ac(0) and to the traditional parameters of Regge theory. 
If we have 

ac{0) + a'^m] = J (36) 
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and J = 2n + 2, then we need 



2a,(0) + 2 = a,(0), 



and 



2a' = a' 



(37) 



In the Regge hniit (where we keep only the leading s behavior), the amplitude from the exchange of the lowest 
particle on the trajectory will simply be 



A 



a'J'hxKt - ml) 

where we have used the fact that the mass of the lowest particle on the trajectory (the spin-2 particle) is 

ac(0) 



2 2 



(38) 



(39) 



If we take the Regge limit of the full amplitude, however, we find 



A 



r 


1 - 


ac(t)' 
2 




r 


a 


c(t) 

2 


- 1 - 


X 



-iizctc (t) /2 



a's ^ 



(40) 



(where we are now using the characteristic parameters ac(0), a'^ of the Regge trajectory). Note again that the factor 
/(ei) contains all information about the incoming and outgoing particles. We can thus relate this amplitude to the 
amplitude in eq. ( 38 1 by replacing the glueball exchange factor f_]^2 with a Reggeized Pomeron propagator 



t — TO^ 





^ 2 




2 r 




-e 



, ^ Qe(t)-2 
-i7rQe(t)/2 f ^ 

2 



(41) 



The factors of /(e^) have indeed cancelled. We can now find the full Pomeron contribution to proton/proton scattering 
by applying this same replacement rule to the graviton propagator in ( 29 ) . The proton/proton differential cross section 
becomes 



da 
~dt 



corresponding to the invariant amplitude 



\^A\t)T'^[^xW^ 


■-, ae(t)' 
^ 2 


7rr2 







2 



2ac(t)-2 



(42) 



A{s,t) = AsX^A\t)e 



-i7rQc(t)/2 



r[- 


-x]r 


1 - 


2 


r 


"Qc(t) 
2 


- 1 


- X 



2 



o(t)-l 



(43) 



This form provides a model for the differential cross section, total cross section and p parameter of either pro- 
ton/proton or proton/anti-proton scattering at very high center of mass energy, where the process is dominated 
by Pomeron exchange. This prediction is relatively model- independent, relying only on the structure of the closed 
string amplitude and the assumption that the graviton (which must be present in any dual theory) couples to the 
energy-momentum tensor. It depends on four parameters: the two trajectory parameters ac(0) and a'^, the glueball- 
proton-proton coupling strength A, and the dipole mass Af^ (for small t we can approximate the form factor A{t) 
with a dipole A{t) = (1 — i/A/J)^^). Of these unknowns, the coupling A, the dipole mass Md, and the glueball mass 
— {2 — q;c(0))/q!^ are all present in the low-energy process involving the exchange of the lowest glueball. That is, 
if we know the low-energy process, which we can compute in the supergravity limit, then the only dependence on the 
full string theory lies in determining the trajectory slope a'^. 

We now present two approaches for fixing the four parameters of the model: (1) we calculate three of them in a 
specific dual model, and (2) we compare these results to least-squares fits of our model to scattering data. 



10 



IV. COMPUTATION OF PARAMETERS IN THE SAKAI-SUGIMOTO MODEL 



After briefly reviewing the Sakai-Sugimoto model |14j . we compute the parameters rUg, A, and Aid- We note first 
that the Sakai-Sugimoto model depends on three quantities, Mkk, Sym, and Is- The string scale Ig does not appear 
in the low-energy supergravity limit, only in the full string theory. The other two parameters are arbitrary a priori, 
but they may be fitted using the p mass and the pion decay constant [T3]. We find excellent agreement between irig, 
Md, and A as computed in the fully fixed Sakai-Sugimoto model, and the values determined by fitting to pp and pp 
scattering data. 



A. Sakai-Sugimoto Model 



The Sakai-Sugimoto model [Hllin!, is a top-down QCD dual: it relies on a brane construction in lOd supergravity 
to produce the salient features of strongly coupled QCD, such as confinement and chiral symmetry breaking. Nf 
flavor I?8-branes are placed in the background generated by Nc Z)4-branes. Closed string (or bulk) excitations are 
dual to QCD glueballs, while open strings living on the probe branes and transforming in the adjoint of a U{Nf) 
symmetry are dual to scalar and vector mesons. 

The DA background is defined by the following metric, dilaton, and Ramond-Ramond three-form C3 (with F4 = 

ds' = (^"' {v^.dx^dx'^ + f{U)dr') + + U'dnl^ , (44) 



*--(^r. /<-)-if 



A radial coordinate U and a unit S"* parametrize the directions transverse to the 1)4- branes. d^l^ is the metric on 
the unit S"*, which has volume form 64 and volume V4 = Stt^/S. = TrgsNJs- The D4-branes are extended in the 
T and x^-directions, for 11 — 0,1, 2, 3. The r-direction is made periodic with t ~ r -I- 2ttAIj^]^. The radial coordinate 
U must now be bounded from below {U > Ukk) to avoid a conical singularity. In order to break the remaining 
supersymmetry, we impose anti-periodic boundary conditions on the fermionic modes so they acquire masses of order 
Mkk- 

It is often useful to work with the scale and the effective 4d coupling constant of the Yang-Mills theory. In terms 
of Ukk and i?, Mkk = 3C/]/^/(2i?3/2), ^ 2TTMKKgsls- 

Placing Nf D8 branes in this background produces flavor degrees of freedom. In the probe limit {Nf ^ Nc) the 
backreaction of the D8-branes with the _D4-brane geometry is negligible. The flavor branes assume a nontrivial profile 
in the {U,t) plane, and are fully extended along the and the directions. 



B. Open and closed string spectra 



Mass spectra of excitations coming from bulk and brane modes can be determined by perturbing the supergravity 
and the brane (DBI) actions, respectively. Computations of the glueball (i.e. closed string) spectrum for this back- 
ground geometry were performed in 'Al\ and |42j . We briefiy review the treatment of [41 using the conventions of 
|14j . and cite relevant results. 

The lOd bulk field content consists of the graviton, Hmn, a dilaton ip, an NS-NS tensor Bmn, and R-R one- and 
three-forms Cm and Cmnl- Neglecting any dependence on the S*^ transverse to the Z)4-branes, and ignoring all but 
the lowest KK modes in the compactified r direction essentially reduces the problem to five dimensions: {x^,U). 
We can classify the states according to their transformation properties under S'0(3) rotations in the physical space 
directions {xi,X2,xz) of the field theory. The graviton in particular gives rise to a 2++ state (/iy), a 1 ^ state 
(/lir), and a O''"*' state (/irr)- The coupling of the bulk fields to the boundary gauge theory and the parity- and 
charge-conjugation invariance of the overall action determine the parity and charge quantum numbers of the 4d field 
theory states. 

Now consider the standard supergravity action 

S=^Jd''x^i^-'^R + A{V4>Y)^^^^Fl^ (45) 
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where kiq is the lOd Newton constant. We introduce perturbations around the background metric by taking Gmn = 
Gmn — fiMNj where Gmn is the background metric from eq. (44l. Varying with respect to hij are spatial Lorcntz 
indices) , we have the equations of motion 



1 /9/ 

2 V 2 



+ Wf ] h,, + (/ + Uf) UdulHj + fU^d^h,, 



u 



(46) 



The prime denotes differentiation with respect to U, and is the 4d momentum of the mode with hij{q,U) = 
J d'^xe~^'^^hij{x,U). We work in a gauge where hoM — 0, dihAB ~ and retain the traceless (spin 2) piece of hjj. 
Writing the lOd perturbation hij{x, [/) as a tower of resonances, 



/i.,(x,c/) = ^/.i;'(x) 



3/2 



TniU) , 



(47) 



the equation of motion becomes an eigenvalue equation for the modes Tn{U) with q^ = rn^, where to„ is the mass of 
the n*'' resonance: 



(48) 



As discussed earlier, we work in a strict Regge limit, only taking into account contributions from the leading Regge 
trajectory, and not from the daughter trajectories. We will therefore use only the lightest mode in the KK tower of 
2++ glueballs. For simplicity, we define T{U) = T'„([/), hij = h\^\ and = ml [63] . The mass of the hghtest 
glueball is proportional to the lowest eigenvalue. 



ml = 1.57 M^^ , 



(49) 



which agrees with the result derived in [IT] . 

Having computed the mass of the lightest spin 2 graviton mode, we must now normalize its wavefunction to yield 
a canonical kinetic term in the effective 4d action. The graviton kinetic term comes from 



S 



1 

2k^q 



which we expand to quadratic order in h^i,: 

1 



Ge-^'I'R^ -V-Ge-^'''G''''&^G"'^ 
2 



1 



(50) 



(51) 



assuming again that the graviton is traceless. Using the expansion in (47 1 and writing the integral in terms of the 
dimensionless ratio U/Ukk, 



357r2 



U 



T^{U)d 



Ukk 



dahsudnh 



1 



J I I'll 



The coefficient of the kinetic term for the 4d spin 2 mode is 



= 



357r2 



U 



T^{U)d 



U 



Ukk J 



3^712 



(52) 



(53) 



where T{U) has dimensions of length and It has dimensions of length squared. Rescaling T{U) T{U)/Nt yields 
a canonically normalized kinetic term. 

Now we turn to the open string spectrum on the Z38-branes, given by the leading DBI action 



S 



DBI 
D8 



= -Tg / d^xe-''' tr V- det {Gmn + 2Tra'FMN) 



(54) 



where Tg = (2Tr)~^lj^ is the brane tension, Gmn is the pull-back of the background metric onto the brane stack, and 
Fmn = dM^N — On Am — *[^m, is the non-abelian field strength of the U{Nf) gauge fields living on the branes. 
We have used the convention ti-{T°-T'') — S"'''/2 for the gauge group generators. 
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By extremizing the DBI action without gauge field fluctuations {Am 
DS-brane stack in the {U, t) plane: 



r(C/)=±C/„V(C/o)^/' 



dU 



0), we can determine the profile of the 



(55) 



where Uq is a constant of integration. The geometry of the branes explicitly realizes chiral symmetry breaking. The 
radial coordinate U corresponds to the energy scale of the dual field theory. Near the UV boundary {U oo), the 
solution exhibits chiral symmetry: it resembles a pair of parallel I?8- and -D8 stacks separated in r by some distance 
L, with 4d modes transforming under [/(iV/)^ x U{Nf)ii. As U (the energy scale) decreases, the D8 and D8 branes 
curve toward each other, until they meet at U = Uq, breaking the chiral symmetry of the two independent UV brane 
stacks to U{Nf)v- Like [2] , we focus on the solution where = Ukk, and the D8 and D8 lie at antipodal points 
on the T circle. It will prove convenient to parametrize the direction along the probe branes in {U,t) plane with 



± 



/ U 



\Uk 



K 



1/2 



(56) 



where Z € (—00,00) such that Z — > ±00 are the left (right) UV boundaries. 

The gauge field fluctuations in the (/i, Z) directions on the D8-brane give rise to towers of 4d vector and axial- vector 
meson states in the field theory. Assuming no dependence on the coordinates, the DBI action (54) becomes 



'r>8 



d'^xdZ tr 



K{Z)-^'^7^^^''7^P-F,.„F, 



MIkK{Z)t^^''F^zF,z 



with 



2167r3 



and K{Z) ^l + Z"^ 



(57) 



(58) 



In order to ensure that the mass and kinetic terms of the Ad action are normalizable, we must have the field 
strengths {F^z, Ffiv) — > as Z ^ ±00. We can choose a gauge where Am itself vanishes at large Z. In order to 
more conveniently realize the (axial)vector meson spectrum and the chiral Lagrangian for the pion modes, we follow 
[H] and transform to a gauge where Az = : 



A^(x^z) 

Az(x'^,Z) ^0 , 
where the gauge transformation g has the form of a Wilson line, 



4'^ = .9^M.9"' + Wd^^g-' 



g{x^,Z)^Pexp 



[ dZ'Azix^'^Z) 
Jo 



(59) 
(60) 



(61) 



A^f^'' now splits naturally into a normalizable piece, gAfj^g^^, which gives rise to the (axial) vector meson spectrum, 
and a nonnormalizable piece igd^g~^. 

Focusing on the non-normalizable modes in the DBI action, we arrive at the action of the 4d Skryme model whose 
solitonic excitations we identify with baryons. The non-normalizable piece of A\^\ igd^g^^, changes the boundary 
conditions on the gauge field such that 



liin 4''H2;'^'±°o) = igdf,g 



^e±(^^)5pC;'(^'') (62) 

where G U{Nf)L and ^_ G U{Nf)fi. Imposing this behavior as a boundary condition, we can express ^If"* as 

00 

= *e+a^e+V^+ + ^i^d,i-yj- + J2 B^.{x)^''Hn{Z) (63) 

n=l 

with 



1 ± — arctan(Z) 

TT 



(64) 
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The nonnormalizable functions ijj±{Z) satisfy the gauge field equations of motion with boundary conditions ■i/j±{Z 
±00) = 1 and ^±{Z Too) = 0. From £,±{x^) we can construct the appropriate chiral field U{x^) = 
transforming as J7 — > under {g+,g-) G U{N f)LxU{N f)^. Following [T3] , we use the residual gauge invariance 

to further fix ^_|_ = U{x^^) and ^_ = 1 so that 

A^^-* — iU^ di^iUip^ + normalizable modes (65) 

We have 5d field strengths 

F^^'J = - [U-\x^^)d^U{x^^),U-\xnd,U{x^^)] ^+{^+ - 1) (66) 
Pzi = iU-\x'')df,U{xndz^P+ . (67) 
The DBI action now explicitly produces the 4d Skyrme model 



5* 



DBI 



K [ d'^xtr^A{U-^{x^')d^U{x^')f + B[U-\x^')^^,U{x^'),U-^{x^')^^U{x^')Y'^ (68) 



with parameters A and B defined by the overlap of the non-normalizable modes with warp factors associated with 
the background metric: 

A = MIj^ J dZ{l + Z^) {dzip+f = 0.318 (69) 
B=\J dZ{l + Z^y^^^Pl (V'+ - 1)^ = 0.078 . (70) 
To match the Skyrme Lagrangian with U{Nf) generators normalized to tr(rT'') = (5°V2, 

Sskyrme = J d^xTv (^ll{U-'d^Uf + ^[U-' d^U, U'' d^?^ (71) 



we identify 



= AkA ^ = 32kB . (72) 



It is easy to show that taking U{x^) — exp[— 2z7r(a:'')T''//Tr] and expanding the action to leading order in the pion 
field indeed gives a canonical kinetic term for the pion. This is the form of the famous Skyrme model |43l 144] . in 
which baryons appear as solitonic configurations of U{x^). The second term in the Lagrangian (the "Skyrme term") 
stabilizes solitons of finite size. Inchiding the u) meson and a gauged WZW term can also be used for this purpose 

SSI- 

We are now in a position to fix the remaining free parameters from experimental data. Following |14] . we use the 
p meson mass rup to fix M^k ~ 949 GeV and the pion decay constant /^r — 93 MeV to fix k = 7.45 x 10"'^ [03]. It 
should be noted that this is the crudest way to set the model parameters, and is intended only to yield a heuristic 
estimate of the values we can derive from the Sakai-Sugimoto model. More accurate results could be obtained by 
fitting to multiple real world parameters (such as several more (axial) vector meson masses). An analysis of this type 
is conducted in the "hard- wall" model of 115 . 



C. Predictions for Free Parameters 



Having detailed the various ingredients of the Sakai-Sugimoto gravity dual, we can now make predictions for three 
of the four parameters appearing in our ansatz for small t, large s proton-proton scattering: the ratio ~ao/a', the 
dipole mass Md in the gravitational form factor, and the coupling A between the proton and the 2++ glucball. The 
first two quantities we take from the existing literature, and calculate the third directly. 



1. Using MxK — 949 MeV and (49 1 we estimate rUg « 1.49 GeV. This value is significantly lower than the 
lattice result, mg-iattice = 2.40 GeV [46 . In the next section, we find that our smaller value of more closely 
approximates the ratio —ao/a' we find by fitting our model to pp and pp-scattering data. 
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2. We model protons as 4d Skyrmions, for which the matrix elements of the energy momentum tensor (EMT) have 
been computed explicitly |38j. The form factors can be related to the components of the static EMT in Breit 
frame (that is with = 0): 

T,Ar,s) = — y —e^^-{p',s'\T,^{0)\p,s) (73) 

where the proton spin polarizations are defined to be equal in the respective protons' rest frames and equal to 
s = (0,s). The form factors are given by 

A{t) + B{t) + j{A'{t) + B\t)) = j d^re-'''-'e'^''s,rjTr,s (74) 

A{t) + [A{t) + 2B{t) + 4C{t)] = - f d3re-''-'^Too(r, s) (75) 
4rn-^ m J 

C'it) + |c'(t) + ^C"(i) = 1 d3re- ''T,,(r) [rW= - ^5'^) (76) 

with m the nucleon mass, and primes denoting differentiation with respect to t. Evaluating the Skyrme EMT on 
the hedgehog solution of |44j . U = exp[jr • fF(r)\ with the radial function F(r) chosen to minimize the soliton 
mass. [55] determine the form factors explicitly in the large Nc limit. A dipole form approximates A{t) well for 
up to \t\ < 0.8 GeV^, with dipole mass Md = 1.17 GeV. This value for Md is in good agreement with the value 
obtained by fitting to data, as presented in the next section. 

A more rigorous analysis would treat the protons as 5d solitons stabilized by vector mesons via the Chern-Simons 
term. This is beyond the scope of the present work and we take the ordinary 4d Skyrme model to be sufficient 
to provide a heuristic estimate for M^- 

3. We now compute the glueball-proton coupling in the Sakai-Sugimoto model, and find that the glueball indeed 
couples primarily to the 4d energy- momentum tensor. 

Let us consider graviton couplings in the DBI action, h^i,. Keeping only couplings linear in h, 



Sd8 = -Ts J d^xe * tr \J -det{GMN - hMN + ^na'FMN) 

= -^Ts{27raT J d^e-^VdC^^'^&^Kpix, z) tr [g-^'F^^F^s + G^^ F^^F, 

= -i^j d'^xKfj{x)tr {Ah{U-^d°'U){U-^d^U) + Bh [U-^d''U,U-^dpU] [U-^d^U,U-^dPU]) + . . . 
= -hapG"^ + ... (77) 



To arrive at the third line, we inserted the expressions for the field strengths in terms of U{x) (eq. 66 1, and 
hap^x, Z), the puUback of the lightest graviton mode onto the branes. The coefficients Afi and B^ are given by 
the overlap integrals 

A„ Jl^il^ r ,ZK(zm^OM^ = 13.42^ (78) 

Bn --nP^ r rf^^(^)-^/'^V4(^+ - 1)' = 7.64 . (79) 

Nc^9ymMkk gvMNp'MKK 



Let us compare the tensor G°'^ in (77 1 to the energy momentum tensor (EMT) of the 4d Skyrme model, 

To^fj = Ktr {A{U-^daU){U-^dpU) + 2B [U^^daU, U^^dpU] [U^^dpU, U^^dfU] ] . (80) 

We can extract the coupling constant A as the ratio of the coefficients of either of the two terms in Ca(3 and 
Taf3 , or some linear combination of both. For example we can consider 

^o/3 ^ ^ ^ fj^^^ _ '^BAjA ^y-iga^^ U-^OpU] [U-^d^U, U-^d^U] (81) 
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where 



Ah 



Or alternatively 



, = 42.18 ^ = 3.61 GeV"^ (82) 

^ MKKgVMNl' 

Bh - - 0.09 GeV-i (83) 



C^P = ^T'^P + K(Ah- ^] {U-^d"U){U-^d^U) (84) 

2l3 \ 2_D J 



where 

Bh .„ „„ 1 



48.86^75 =4.19 GeV"^ (85) 

^h~^)= -0-17 GeV . (86) 



2B 



(87) 



Assuming the relative contributions of the two terms (kinetic and Skyrme) are of the same order, the deviation 
of Cq/3 from the EMT amounts to a few percent of the value of A. We therefore estimate A ~ 3.9 ± 0.3 GeV^^, 
in excellent agreement with the fits discussed in the next section. It should be kept in mind however that the 
Skyrme model itself is only accurate to within ~ 20% so the value of A derived from the Skyrme model may 
deviate from its actual value by somewhat more than this estimate. 



V. DATA FITTING AND COMPARISON 



A. Regime of validity 



Before using our model to fit experimental scattering data, we briefly discuss the data we use, and limitations in 
applicability of our model as a function of s and t. We fit the differential cross section da/dt for a variety of values 
of s and t (shown in Fig. |2|. There are many complicating factors which limit the validity of our model, the most 
obvious being that we have worked in the strict Regge limit, neglecting corrections suppressed by powers oi t/s. For 
all of the data we use, t/s < .001. Other effects, such as Coulomb contributions, perturbative QCD effects, and the 
contribution of Reggeon trajectories are not so easily discarded. We discuss each of these in some detail below. 



1. Coulomb Contributions 



For very small values of |t| (regardless of s) the Coulomb interaction makes a significant contribution to the 
amplitude. This contribution will be largest at \t\ ~ 0.002 GeV^ and is negligible (for our purposes) by t ~ 0.01 GeV^ 
[UliH]. We will use |i| = 0.01 GeV^ as a lower cutoff in t. 



2. Lower Regge trajectories 



For large values of s, the total cross sections converge because the exchange of a Pomeron does not distinguish 
between particles and antiparticles. For small values of s, however, there are contributions to both proton-proton 
and proton-antiproton scattering from other Regge trajectories. In particular, the lower Regge trajectory is actually 
a pair of exchange degenerate trajectories, one consisting of even spin particles and the other of odd spin particles. 
For proton-antiproton scattering the Reggeon contribution is larger because these two trajectories add, whereas for 
proton-proton scattering they work to cancel each other out. We can estimate how large the contributions from the 
next Regge trajectory will be by looking at the total cross sections for proton-proton and proton-antiproton scattering. 

We fit the total cross sections with the functions 



cTtotipp ^ pp) = PsP + Qs'^ 



(88) 
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FIG. 2: Data from proton-proton and proton-antiproton scattering, with a dot for every data point in (logs, log space. A 
given experiment is generally run at fixed s, for a range of values of t, so it creates a vertical series of data points. Below the 
diagonal line t/s corrections make up less than 0.1%. The horizontal line near the bottom lies at \t\ = 0.01 GeV^, which we are 
using as a lower cutoff. Below this line there are significant effects from Coulomb scattering. The horizontal line near the top is 
at \t\ = 0.6 GeV^ which we use as an upper cutoff. Along the log s axis we mark the values at which the Reggeon contribution 
to proton/proton or proton/anti-proton scattering is 1%, 5%, and 10%. 




FIG. 3; Log-log plot of total cross sections (in mb) for proton-proton and proton-antiproton scattering as a function of s (in 
GeV^). Note that the cross sections converge as s grows large. 
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^tot {pp ^ PP) = PsP + Rs" (89) 

and find best fit values 

p=.08, g=-.46, P = 21.3, Q = 53.2, i? = 103.6. (90) 

We can use the optical theorem to relate the total cross section to the differential cross section at t = 0, and use 
this in turn to estimate the size of the first contribution from the lower trajectory. 

^{ppipp)) - P's'^ (l + ^.^-^ + ...). (91) 

Based on this functional form, the magnitude of Reggeon contamination in proton/anti-proton scattering at ^/s = 31 
GeV is about 22%, while at = 1800 GeV it is about 0.3%. In FIG. [2] the marks on the logs axis show where 
Reggeon contamination to proton-antiproton and proton-proton scattering is 1%, 5% and 10%. For the lower center- 
of-mass energy data, the effect is fairly large. We could account for this by adding to our model a term corresponding 
to Reggeon exchange. For the present treatment, however, we simply add the amount of Reggeon contamination for 
a given value of s to the experimental error associated with each data point, thus weighting our fit towards the higher 
energy data. 



3. Perturbative QCD and the hard Pomeron 

For sufficiently small values of \t\, it is reasonable to think of the Pomeron as a trajectory of confined states 
(glueballs). This is the "soft Pomeron" on which our model is based. As \t\ increases, however, we eventually enter 
the regime of perturbative QCD, outside our model's regime of validity. It is not clear on theoretical grounds exactly 
where this transition occurs; we instead attempt to determine its location empirically, by examining the data. Let us 
assume that the Pomeron contribution to the differential cross section is of the form 

— = F(t)s2"^(*)-2 (92) 
at 

for all values of t, where F(t) is some unknown function. Differential cross section data is typically analyzed for a 
range of t values at a fixed value of s. Suppose instead we consider a fixed value of t for a range of values of s. The 
trajectory adt) is therefore the slope of log{da/dt) plotted as a function of log{s): 

log ^ = log F{t) + {2a,{t) - 2) log s. (93) 

By fitting for the slope of this line at a range of values of t, we can get a reasonable picture of the function adt). 
Referring again to FIG. |2] we can see that taking sets of the data at fixed values of \t\ is difficult with the extant 
data, and generally these sets will only have between 3 and 5 data points each. This grouping of the data would 
not yield reliable statistics, so we do not use it to fit adt) directly, but consider it a reasonable estimate of where 
the transition occurs between the regime where the soft Pomeron accurately characterizes the exchanged degrees of 
freedom, and the regime where it does not. A graph of adt) as a function of \t\ is shown in FIG. |4] We can see 
that for < |t| < 0.6 GeV^ the trajectory matches what we expect for the soft Pomeron. There is a clear transition 
at |t| = 0.6 GeV^, where the slope suddenly becomes much less steep; above \t\ > 0.9 GeV^ the behaviour is clearly 
nonlinear. We therefore attach no particular significance to the shape of the plot in this region, choosing instead to 
impose an upper bound of 0.6 GeV^ on the values of \t\ we consider in our fits. 



B. Comparison to a "Photo-Pomeron" model 



In order to provide a frame of reference for the success of the model described above compared to the existing 
literature, we briefly review a commonly-used model for single Pomeron exchange due to Donnachie, Jaroszkiewicz, 
Landshoff, Polkinghorne and others ([4^], [5U], [5T]). This is only one of many possible models of various degrees of 
complication, but serves as an illustrative example because it has the same structure as our model, with the difference 
that it relies on the electromagnetic proton form factor rather than the gravitational form factor. Other models based 
on different assumptions include the impact picture model of [55], the multi-component model of the Durham group 
(see [53] [SJ for a recent discussion) and the eikonal model of Block, Halzen and collaborators (see e.g. [S5] and 
references therein). 
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FIG. 4: The function ac{t). In the region < |t| < 0.6 GeV^ the trajectory matches that of the soft Pomeron. Outside this 
region we assume that perturbative QCD effects make our treatment inappUcable. 



The electromagnetic- type Pomeron coupling (reviewed in \52\ ) draws inspiration from the additive quark rule: the 
(experimental) fact that the ratios of total cross section equal the ratios of the numbers of valence quarks present in 
the scattered hadrons. Positing that the Pomeron couples to constituent quarks individually as 7'^ reproduces this 
observation. The form of the Pomeron-quark coupling is then assumed to be identical to the photon coupling except 
that the Pomeron is C = +1 rather than C — —1. The form factor Fi{t) involved in the exchange is assumed to 
be identical to the electromagnetic form factor. For \t\ < 1 GeV^, a dipole approximation to Fi{t) (from electron 
scattering data) gives 

Ami - 2.7m 1 

Using this single-Pomeron exchange model, the unpolarized pp (or pp) cross section becomes 

da ^ (3/3pFi(t))4 / . y-"^*)-^ 

dt An \sq) ' ^ ' 

with So ~ 1 GeV'^ the characteristic scale of the problem. The key difference between this model and ours lies in the 
form factor: gravitational in our case, electromagnetic in the case of the "photo-Pomeron." 



C. Fits to Scattering Data 

We perform standard least squares fits |65j to the differential cross section data for pp and pp scattering using the 



form (42 1, allowing the parameters aoj Q^', M^, and A to vary. We also perform an identical fit for the photo-Pomeron 
model, allowing /3p, and a' to vary. The data is taken from the Durham HEP database (http://durpdg.dur. ac.uk| 
with 0.01 < t < 0.6 GeV^ and 30.4 < y/s < 1800 GeV. As our model does not take into account the effects of Reggeon 
exchange, we simply estimate the contribution of Reggeons to da /dt at particular values of i/i, and add the result in 
quadrature to the experimental errors. There is significant disagreement between two data sets at -^/s — 1800 GeV 
produced by the CDF and E710 experiments, as shown in FIG.jS]. Rather than choosing one or the other set explicitly, 
we perform all fits using both, just CDF, and just E710, with results as displayed in Table |T] 

Our model clearly produces the smallest value of /d.o.f. when only the CDF data set is included, but fares better 
than the photo-Pomeron model in all cases. FIG. |6] shows fits of our model to the differential cross section data. 

We can also compare our predictions for the total cross section and p = YleA{t = Q)/lmA{t = 0) (based on the 



best-fit parameters determined from the differential cross section) to data. Applying the optical theorem to eq. ( 40 1 
yields the total cross section 
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FIG. 5: A log-linear plot of the differential cross section at -^s — 1800 GeV from the CDF and E710 collaborations, as a 
function of t (in GeV^) . 



Gravitationally coupled Pomeron 



both data sets 


just E710 


just CDF 


ac(0) = 1.076 ±.0016 
a'^ = .290 ± .006 GeV"^ 
Md = .983±.016GeV 
A = 4.28 ± .03 GeV-' 

d.o.f 


ac(0) = 1.074 ±.0016 
a'e = .286 ± .006 GeV"^ 
Md = .970 ± .016 GeV 
A = 4.31 ± .03 GeV"' 

ttS = 1-41 


ac(0) = 1.086 ±.0016 

= .300 ± .006 GeV"^ 
Md = 1.02 ± .016 GeV 
A = 4.14 ± .03 GeV"' 

y, = 1.26 

d.o.j. 


Electromagnetically coupled Pomeron 


both data sets 


just E710 


just CDF 


ac(0) = 1.076 ±.0013 

= .289 ± .003 GeV"^ 
f3 = 1.858 ± .016 GeV"' 
= 1-97 

d.o.f 


ac(0) = 1.075 ±.0013 

= -289 ± .003 GeV"^ 
P = 1.877 ± .016 GeV"' 
= 1.66 

d.o..f 


Qc(0) = 1.082 ±.0018 

= .289 ± .003 GeV"^ 
/3 = 1.801 ± .020 GeV"' 
= 1.79 



TABLE I: Fits to differential cross section data including both E710 and CDF ^/s — 1800 data sets, just E710, or just CDF. 



where in terms of the best fit parameters (using only CDF data) we find C — 21.325 and ao — 1 = 0.085. Performing 
an expHcit fit to total cross section data, we find Cfn = 21.097 and (ac(0) — 1) fu — 0.086, in excellent agreement 
with the computation. 

Because we neglect Reggeons, (40 1 predicts a constant value for p as p = — cot — tt) = 0.136 (where again we 
use the values of the best-fit parameters from the differential cross section) . This value agrees well with the data at 
large ^/s (see Fig. [7l, where the Reggeon contribution is minimal. 

We have shown that our mechanism for Pomeron exchange in pp and pp scattering fits experimental data quite 
well. The best fit parameters from fitting the data also compare favorably with our estimates for ac(0)/a^. A, and 
Md computed in holographic QCD. The Sakai-Sugimoto model predicts the mass of the lightest spin 2 glueball 
irig = 1.485 GeV, while the value produced by fitting the form of the differential cross section to scattering data 
yields nig-fit — 1.745 ± 0.035 GeV, is within 15% of the computed mass, though not with the statistical error bars 
determined by the fit. The gravitational dipole mass Mj, computed in the Skyrme model has value = 1.17 GeV, 
which deviates from the fitted value M^-fit = 1-02 ± 0.016 by about 15% as well. As the Skyrme model predicts 
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.300 ± .006 GeV"^ 


Md = 


1.02 ± .016 GeV 


A = 


4.14 ± .04 GeV"^ 




1.26 
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FIG. 6: A log-linear plot of the best fit to the scattering data using just the CDF data at ^/s = 1800 GeV. The differential 
cross section is in mb/GeV^ and t is in GeV^. Note that at lower center of mass energies the fit is less successful, most likely 
due to a greater contribution from Reggeon exchange. 
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FIG. 7: Data for p as a function of logs compared to the predicted value. 



masses only to an accuracy of about 20%, the fitted value lies within the expected uncertainty of the computed dipole 
mass. Finally, the coupling constant computed from holography to be A « 3.90 ± 0.3 GeV~^ agrees extremely well 
with the best fit value of 4.14 ± 0.04 GeV"^ 

We should note that the values of per d.o.f. we obtain are not as good as those of fits cited in the current PDG 
which typically fit to a leading log^ s behavior and use more sophisticated filtering of the data than we have done 
[521 [SHIP]- 

Finally, we can use our model for pp scattering to make a prediction for the differential and total cross sections to 
be observed at the LHC (at ^/s — 14 TeV). The differential cross section da/dt is plotted in Fig. [sjas a function of 
the momentum transfer, t. We predict a total cross section of atot = 109 ± 4mb. 
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-6L 

FIG. 8: A log-linear plot of the predicted LHC differential cross section (in mb/GeV^) as a function of t (in GeV^) at 
y'i = 14 TeV. We have used the best fit parameters from the data set including only the CDF ^/s — 1800 GeV data. The 
range shown is generated by the errors in the fit parameters. 

VI. CONCLUSIONS AND FUTURE DIRECTIONS 

We have used AdS/QCD to construct a general model for Pomeron exchange in Regge-limit p—p and p—p scattering. 
In order to do so, we assumed that string scattering ampUtudes have the same structure in weakly curved space as they 
do in flat space, but that the values of the Regge trajectory parameters and the masses of excitations are modifled. 
This implies that the curved space Regge regime amplitude factorizes into a piece that characterizes the interaction of 
the scattered particles with the exchanged trajectory, and a piece that corresponds simply to the exchange of a closed 
string (the Reggeized propagator). Furthermore, the coupling of the external states to the exchanged trajectory is 
described entirely by their coupling to the lowest mode (the 2++ glueball). 

Using these principles, and generic properties of QCD duals, we were able to identify the form factor A{t) at the 
Pomeron-proton vertex as coming from the matrix element of the energy-momentum tensor, with the strength of the 
coupling given by an overlap of graviton and pion wave functions. In a particular hQCD dual, the Sakai-Sugimoto 
model, we computed the coupling directly and showed that it agrees with one determined by fits to experimental data. 
We also developed a method for extending the amplitude for exchanging the 2++ glueball to include the exchange of 
the entire glueball trajectory. 

Though our treatment offers several advantages compared to previous approaches to this problem, there are some 
ways in which our methodology could be improved or extended. At the level of numerical analysis, our errors did not 
take into account systematic errors that might drive an entire data set at a particular -y/s either up or down. More 
sophisticated fitting techniques would also filter outliers out of the data, which could significantly improve the we 
find. 

On a theoretical level, one could certainly extend the regime of validity of our model to lower values of ^/s by 
modeling Reggeon (open string) exchange as well as Pomeron exchange. Our treatment of the proton in the dual 
model was also rather simplistic. Computing the proton form factors and proton-proton-glueball coupling using the 
duals of baryons as 5d instantons in the Sakai-Sugimoto model [501 EI] rather than via the simple 4d Skyrme model, 
would yield a more accurate holographic picture; one should also take into account that the 5d solitons may be 
stabilized using vector meson modes. As mentioned above, [35] recently presented a treatment of protons as fermionic 
fields in bottom-up holographic models, which included results for the electromagnetic and gravitational form factors. 
It would be interesting to incorporate these techniques, possibly coupled with the Skyrmion picture of protons, into 
our model. 

We used the four-point sphere amplitude to model the "Reggeized propagator." In the Regge limit this amplitude 
indeed consists of the t-channel exchange of a closed string, to which the Pomeron is dual. However, it is not clear that 
the incoming and outgoing particles (the protons) are themselves dual to closed strings. As they may be considered 
as either to solitonic configurations of open string modes (pions) or as wrapped D4-branes, we might need to consider 
a more complicated amplitude to accurately refiect the structure of the scattering process. This is a difficult problem, 
but additional insight might be gained from tt — tt scattering, where the string dual should be an annulus amplitude. 

Finally, there is the difficult issue of corrections for large s. The Froissart bound indicates that at some high s, the 
behavior s°'('^^'^ for the total cross section must be replaced by a function that grows at most like log^ s. There are two 
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possible sources for large s corrections to our model: string loop corrections, a.k.a. Regge cuts or multiple Pomeron 
exchange, and corrections to the string amplitude from the curvature of the AdS space. It is certainly possible that 
such effects already play a role at energies we consider. More thoroughly examining the effects of spacetime curvature 
in particular would improve the accuracy of our predictions, and would hopefully serve as evidence for the existence 
of a curved-space string dual to QCD. 
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